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ABSTRACT 

We present spectroscopy data for a sample of 14 FUors and 2 TTauri stars 
observed with the Spitzer Space Telescope or with the Infrared Space 
Observatory (ISO). Based on the appearance of the 10 /xm sihcate feature we 
define 2 categories of FUors. Objects showing the sihcate feature in absorption 
(Category 1) are still embedded in a dusty and icy envelope. The shape of the 
10 /im silicate absorption bands is compared to typical dust compositions of the 
interstellar medium and found to be in general agreement. Only one object 
(RNO IB) appears to be too rich in amorphous pyroxene dust, but a superposed 
emission feature can explain the observed shape. We derive optical depths and 



-'^Based on observations with ISO, an ESA project with instruments funded by ESA Member States 
(especially the PI countries: France, Germany, the Netherlands and the United Kingdom) and with the 
participation of ISAS and NASA. This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
a contract with NASA. 
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extinction values from the silicate band and additional ice bands at 6.0, 6.8 and 
15.2 yum. In particular the analysis of the CO2 ice band at 15.2 /xm allows us 
to search for evidence for ice processing and constrains whether the absorbing 
material is physically linked to the central object or in the foreground. For 
objects showing the silicate feature in emission (Category 2), we argue that the 
emission comes from the surface layer of accretion disks. Analyzing the dust 
composition reveals that significant grain growth has already taken place within 
the accretion disks, but no clear indications for crystallization are present. We 
discuss how these observational results can be explained in the picture of a young, 
and highly active accretion disk. Finally, a framework is proposed as to how the 
two categories of FUors can be understood in a general paradigm of the evolution 
of young, low-mass stars. Only one object (Parsamian 21) shows PAH emission 
features. Their shapes, however, are often seen toward evolved stars and we 
question the object's status as a FUor and discuss other possible classifications. 
Two spectra (RNO IB and L1551 IRS 5) show [Fe II] emission lines which are 
attributed to hot and dense material located close to the root of an outflow 
(L1551 IRS 5) or to shocked material directly within an outflow (RNO IB). 

Subject headings: accretion, accretion disks — circumstellar matter — stars: 
formation — stars: pre-main sequence — dust, extinction — stars: individual 
(FU Ori) 



1. Introduction 

By means of mid-infrared (MIR) spectroscopy gaseous and solid state features have 
been observed and analyzed in a variety of astronomical environments. While ground-based 
observations are restricted to certain atmospheric windows (e.g., the N-band around 10 /im), 
spectrographs onboard space-borne telescopes, such as ISO and Spitzer, enable us to study 
a broad wavelength range from the near-infrared (NIR) to well beyond 30 /im. Those in- 
struments fostered and revolutionized our understanding in numerous astronomical fields of 
research. In particular the star-formation community benefitted from studies based on data 
from the space telescopes. The composition of dust grains and their evolution (e.g., grain 



intermediate mass Herbig Ae/Be stars (HAeBes) (e. 


ft' : 


Bouwman et al. 2001: Meeus et al. 


2001: Acke & van den Anckei 


2004 


; \ 


^an Boekel et al. 


2005r), but also in the voune. less mas- 


sive TTauri stars (e.g.. 


Forrest et al. 


2004: 


Kessler-Silacci et al. 2006: Sargent et al. 2006). 



In addition, the ice and dust features of younger and more deeply embedded objects were 
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[e.g., 



Knez et al 



2005 



studi ed fe.g.. IWatson et al.ll2004r) and the ice inventory of molecular clouds was investigated 



Bergin et al.ll2005r ). 



One special sub-group of young low- mass objects are FU Orionis objects (FUors), named 
after the prototype FU Ori. For most objects of this group evidence was found for a tremen- 
dous outburst in optical or NIR light over short timescales (months to years) followed by 
a decline in luminosity typically over several decades. Other objects were included to this 
group as they shared peculiar spectral features, e.g, CO bandhead absorption in the NIR and 
a changing spectral type with wavelength. Most observational data can be well explained 
with highly active accretion disks surrounding these objects, possibly fed with fresh material 
from a remnant envelope (for a review, see, iHartmann fc KenyonI Il996l ). Observations in 
the MIR are thus ideally suited to probe the dusty component of the circumstellar material 
of these objects, either contained in the accretion disks or, additionally, in the envelopes. 
Recent findings from NIR and MIR interferometers support the presence of accretion disk s 
several tens of AU in size (lOuanz et allbood : iMalbet et al1l2005l : lMillan-Gabet et allboOGh . 
However, while the SED and the NIR and MIR visibilities of FU Ori itself can be suffi- 



ciently explained by a simple accretion disk model (iQuanz et al.ll2006l ). iMillan-Gabet et al. 
( I2OO6I ) found that accretion disks alone can not reproduce the SED and observed low K- 
band visibilities of V1057 Cyg, V1515 Cyg and Z CMa-SE simultaneously. They concluded 
that additional uncorrected flux is required possibly arising due to scattering by large dusty 
envelopes. 



While iLorenzetti et al.l (120001 ) presented far-infrared spectroscopy data for 6 FUors ob- 
served with ISO-LW S, a dedicated MIR study of a larger sample of FUors is still missing. 
Larsson et al.l (I2OOOI ) showed the ISO-SWS spectra for 6 FUors for comparison, but no anal- 
ysis was carried out. The ISO data for Z CMa wer e pre sented in publica tions related to 
Herbig Ae/Be stars by lAcke fc van den Anckerl (l2004l ) and lElia et al.l f|2004j). as Z CM a pre- 
sumably is a binary system consisting of a Herbig star and a FUor. IWhite et al.l (120001 ) used 
the ISO observation s of L1551 IRS 5 as input for a radiative transfer model of this source. 
Hanner et al.l (119981 ) discussed ground-based 8-13/xm spectra for four FUors (FU Ori, V1515 
Cyg, V1057 Cyg, V1735 Cyg) and fltted a simple dust model to the data to check whether 



silicate parti c. 



es fro m the interstellar medium (ISM) can reproduce the observed features. 



Schiitz et al.l (120051 ) pubhshed additional ground-based data for another four objects in the 



same wavelength regime. As three objects (Z CMa, V346 Nor, V883 Ori) showed the lO/xm 
feature in absorption they inferred the optical depth from fltting an ISM dust model to the 
spectra. The f o urth object (Bran 76, alias BBW 76) was not analyzed in greater detail. 
Polomski et al.l (120051 ) presented data on RNO IB, Z CMa, and Parsamian 21 and derived 
dust temperatures and optical depths. 



-4- 



To our knowledge , the first MIR spectra of FUors observed witli Spitzer were presented 



m I Green et al.l (120061 ) . Tlie main focus of this publication was an accurate SED modeling 
of FU Ori, V1515 Cyg, and V1057 Cyg. Spectral solid state features were no t analyzed in 
greate r detail. The spectrum of V346 Nor was presented for comparison. In iQuanz et al. 



(120061 ) the first detailed dust composition modeling for the FU Ori spectrum was presented 
and evidence for grain gro wth in the a c cretio n disk was found. The Spitzer spectra of 
RNO IB/IC were shown in lQuanz et al.l (120071 ). However, the dust composition was not yet 
analyzed in detail. 

In this paper we compile MIR spectra for 14 FU Orionis objects observed with ISO 
and/or Spitzer. As up to now only ~20 FUors or FUor candidates are known, this is the 
largest sample of these objects analyzed so far in a single MIR study. Part of the data have 
not been published before. For objects where the spectra show a sufficient signal-to-noise 
ratio, the dust and ice composition of the circumstellar material is investigated. 



2. Observations and Data Reduction 

The mid-infrared spectra we present in this paper are compiled from the archives of 
the ISO§ and the Spitzer§ satellite. The Short Wavelength Spectrograph (SWS) onboard 
ISO consisted of two nearly independent grating spectrometers with a spectral resolution 
between 1000 and 2500 (depending on the band and order) and covering a wavelength range 
from 2.4- 12.0 /im and from 12.0- 45.0 /im, respectively. While the field-of-view (FOV) for 
the shorter wavelength regime was 14" x 20", the FOV for the longer wavelength range 
was 14" X 27", with the exception of the wavelength range between 28.9- 45.0 /zm which 
had a FOV of 20" x 33". With ISO/SWS seven objects classified as FUors were observed 
between April 1996 and October 1997. For one object (00 Ser) data were taken at five 
different epochs documenting a decay in luminosity over a few months. Table [1] summarizes 
the ISO-observations with object names, coordinates, ISO-SWS observing mode and scan 
speed, integration time on target, possible pointing offsets (see below), and the date of 
the observation. For the data reduction Highly Processed Data Products (HPDP) or SWS 
Auto Analysis Results (AAR) were downloaded from the ISO archive for speed 1 and 2 or 
speed 3 and 4 observations, respectively. With the OSIA software package (version 4.0)^ the 
following reduction steps were carried out: For each object the spectra from the ISO-SWS 



-^The ISO archive can be accessed via lhttp://www. iso.vilspa.esa.es/ida/ 
^http: / /ssc. spitzer. caltech.edu/archanaly/archive. html 
ht t p : / / SWS . ster . kuleuven .ac.be/osia/ 
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up- and down-scan were flat-fielded and rebinned. After sigma clipping, the speed 3 and 
speed 4 spectra were de-fringed. This procedure was not required for the speed 1 and speed 
2 data as the HPDP are already de-fringed. Finally, the spectra from the up- and down-scan 
were combined and rebinned to a spectral resolution of 100. In case the resulting spectrum 
showed signs of a pointing offset (e.g., aperture jumps) a correction based on the measured 
beam profiles along the different axes was applied to the raw data and the data reduction 
was repeated. The applied offsets are listed in Table [H 

The Spitzer observations are summarized in Table El The IRS onboard Spitzer offers 
a short wavelength, low resolution module (SL) covering the wavelength range between 5.2- 
14.5 //m, and a short wavelength, high resolution module (SH) going from 9.9- 19.6 /im. The 
corresponding slit sizes are ~ 3.6" x 136.0" (including both SL orders and an additional bonus 
segment connecting both orders) and 4.7" x 11.3", respectively. For the longer wavelength 
part a long wavelength, low resolution module (LL) ranging from 14.0- 38.0 /xm and a long 
wavelength, high resolution model (LH) covering the regime from 18.7- 37.2 /im are available. 
The slits sizes are 10.5" x 360.0" (including both LL orders and an additional bonus segment 
connecting both orders) and 11.1" x 22.3", respectively. Both low resolution modules (SL 
and LL) offer a spectral resolution between 64- 128 (depending on the wavelength) while the 
high resolution modules both have ~ 600. All objects listed in Table [2] were observed with 
the full wavelength coverage, either with a combination of SL-I-LL or with SL-t-SH+LH. 
Two objects (HL Tau and XZ Tau) are not classified as FUors but were part of a small 
Spitzer/IRS map including L1551 IRS 5, and data for all three objects could be downloaded 
simultaneously. Inter estingly, XZ Tau i s a binary system that recently was found to show 



EXor-type variations (jCoffey et al.ll2004l ). i.e., another type of short term eruptions of young 
stars. Thus, a comparison to the FUor data presented here is reasonable. In addition, the 
data of HL Tau enable us to compare the FUor spectra to that of a well-studie d young star 



with a remnant envelope and a highly inclined accretion disk seen almost edge-o n (IClose et al. 



9971) . Part of the Spitzer spectrum of HL Tau was already published in iBergin et al. 
(I2OO5I ). The object V1647 Ori was observed three times within a period of roughly 5 months 
to monitor its brightness as it underwent an eruption beginning of 2004. However, to our 
knowledge, thus far no spectrum was published. 



The data reduction process of the Spitzer data is described in detail in iQuanz et al. 



(120071 ). However, we should mention that the error bar for each individual spectral point 
represents the formal standard deviation from the mean flux of at least two independent 
measurements (two telescope nod position and possibly several observation cycles). Also 
taken into account is the formal erro r of the spectral respon se function. For details on the 



method we refer the reader also to iBouwman et al.l (120061 ). We estimate a relative flux 



calibration error across a spectral order of ~ 5% and an absolute calibration error between 
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orders or modules of ~ 10%. In particular, for the objects RNO IB and RNO IC apparent 
flux density offsets between the SL and the SH part of the spectra, as well as between 
th e short and l o ng w avelengths part of the high-resolution spectra, are already discussed 
in iQuanz et al.l (120071 ) . For the other objects presented here, the discrepancies in the flux 
densities between the SL and the SH part of the spectrum were < 10% and we matched the 
longer wavelength part to the shorter regime by multiplying a scalar factor. O nly Bran 76 



(also known as BBW 76) showed a larger offset of ~ 15% as already mentioned by lGreen et al. 



( 120061 ). For Parsamian 21 (HBC 687) we do not show the SH spectrum between 14 and 
20 fim as the slit of the spectrograph was not centered on the source and significant flux loss 
occurred for which we could not correct. An additional correction between the SH and the 
LH module of the spectrum was required for XZ Tau, HL Tau and L1551 IRS 5 where the 
flux densities of the LH part had to be scaled down by 10- 15%. This offset can be explained 
by the larger aperture of this module which possibly probed additional large scale emission 
from the surroundings of these objects. 



Results 



3.1. General overview 



Figures [HIH] show the complete sample of spectra. For Bran 76 and the third observation 
of V1647 Ori only low-resolution Spitzer data were available. To increase the signal to noise, 
all Spitzer LH data and the SH data of V1735 Cyg were smoothed by a factor of three. 
Most of the yet remaining spikes in these spectra are not real but rather flux jumps between 
the different orders of the spectrographs. As mentioned in the intro d uction of this paper, 
parts of the data shown here were already published: iGreen et al.l (120061 ) presented the 
Spitzer data for FU Ori, V1515 Cyg, V1057 Cyg, Bran 76, and V346 Nor a nd used disk- 



envelo p e models to explain the SEDs. The spectrum of HL Tau was shown by iBergin et al. 



( 120051 ). ISchiitz et al.l ((20051) used the ISO data for Z C Ma and V346 Nor to compare with 
their ground-based data. Finally, iLarsson et al.l (l2000l ) showed the ISO/SWS SEDs of the 
outbursting object 00 Ser and used the ISO/SWS data for RNOIB, Z CMa, V1057 Cyg, 
and VI 735 Cyg for comparison. 

Unfortunately, for most objects the quality of the ISO/SWS data is significantly worse 
than that of the Spitzer observations. Even after the data have been rebinned to a spectral 
resolution of 100, artefacts, i.e. potential emission features that were only detected either 
during the up- or the down-scan, remain in the spectra. In Figure M we show data for 
four objects that were observed with both Spitzer and ISO. In particular, for the objects 
with lower flux levels the noise in the SWS data is significant. The reason for this is the 
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short integration time for most objects which is reflected in the speed parameter in Table [TJ 
In consequence, the ISO data are mainly used for qualitative statements rather then for 
quantitative analyses throughout the rest of the paper. Only the ISO data for Z CMa (with 
a high flux level) and Reipurth 50 (with long integration time) will be examined in more 
detail in one of the consecutive sections. The data on Reipurth 50 is published for the first 
time. 

From Figures [T]1H] it becomes clear that for all objects the flux densities increase toward 
longer wavelengths indicative of warm dusty material surrounding all objects. However, Fig- 
ures m [21 [5l and [6] show that there are striking differences within the group of FUors: While 
some objects show a silicate emission feature in the 10 /xm region, other objects show deep 
absorption profiles. Since other spectral features do further support such a differentiation 
we will in the following distinguish the objects via the behavior of their 10 fim feature and 
discuss the two categories separately in the following subsections. A complete overview of 
the most prominent spectral featu res between 3 and 16 fim is given in Table [31 If we count in 



V883 Ori, which was observed by ISchiitz et al.l (l2005l ). and disregard XZ Tau and HL Tau 



then 9 FUors show the silicate feature in absorption while 6 FUors show silicate emission. 



3.2. Objects with 10 jj,in emission 
3.2.1. Qualitative analysis of the 10 fim region 



In Figure [TOl we compiled all objects showing signs of silicate emission in the 10 fim 
band. To subtract the underlying continuum we fitted a polynomial of first order to the 
flux at 8 and 13 yum0. As we are thus far only interested in a qualitative comparison among 
the various objects, the exact shape of the continuum is not important and higher oder 
polynomial fits did not alter the results relative to each other. For comparison we overplot 
the silicate emission feature of typical interstellar medium (ISM) dust grains, scaled to the 
observed spectra (red, dashed lines in Figure [TOl [Kemper et al.ll2004l ). We note that the 
spectrum of each object is shown twice: While the first plot shows the observed spectrum 
of the source, the second one shows the dereddened spectrum. For the dereddening we 
compiled optical extinction val ues from the l iterature (see caption of Figure [TU]) and used 
the extinction law described in [Mathid (Il990l ) to derive corresponding extinction values for 
the MIR. It is assumed that all of the extinction is caused by material in the line of sight 



''For Parsamian 21 we had to shift the left point of the fit to 6.6 /iin due to the special shape of the 
spectrum (see below). 
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toward the objects and that no self-shadowing effect (e.g., by an inchned accretion disk with 
a certain flaring angle) is present. The continuum was fitted for the dereddened spectra 
separately. It shows that for most objects the extinction is not negligible and that it can 
have substantial influence on the shape of the silicate feature (see, e.g., V1057 Cyg and 
V1647 Ori). For a fair comparison it thus seems reasonable to evaluate the dereddened and 
not the observed spectra. 

At first glance, when comparing the spectra, it is notew orthy that al l emiss ion features 
differ from the typical shape of the typical ISM dust feature. I Green et al.l (120061 ) stated that 
the emission peak for the first four objects in Figure [TO] was close to 9.5 fim and that the dust 
features appeared to be pristin^. We f ind that all features peak longward of 9.7 /xm (the 
typical peak position for ISM type dust; iKemper et al.ll2004l ) and that they show additional 
flux excess compared to the ISM feature at even longer wavelengths. This indicates that 
dust grain processing has already set in. Furthermore, it should be noted that the spectrum 
of Parsamian 21 looks significantly different compared to the others. The most prominent 
characteristic are strong emission bands aroun d 8.2 /xm and probably also around 11.3 /im 
from polycyclic aromatic hydrocarbons (PAHs). iPolomski et al.l (l2005l ) already suspected the 
existence of PAH emission in the spectrum of this source based on ground-based observations, 
but a firm confirmation was thus far lacking. We discuss the spectrum of Parsamian 21 in 
more detail below. 

For a better comparison of the emission features, we plo tted in Figure [TTj t he no rmal- 
ized fluxes of the objects between 8 and 13 /im. 
normalization was done via 



Following Ivan Boekel et al. 



([20051), the 



.(A) = l + 



(A)-Fcont(A) 



> 



(1) 



where Fobs(A) is the observed flux, Fcont(A) is the continuum flux and < Fcont > denotes 
the mean value of the underlying continuum in the considered wavelength regime. This 
normalization ensures that the shape of the emission feature is preserved. From Figure [11] 
it becomes clear that the emission features of Bran 76, FU Ori, and V1515 Cyg are quite 
similar in terms of shape and strength. The features of V1057 Cyg and XZ Tau are far 
less pronounced and much broadeil^. The emission profiles of V1647 Ori and Parsamian 21 
are slightly stronger than those of the other objects and, as mentioned above, the latter 



5 



It should be mentioned that 



6- 



Green et al. ( 20061 ) determined the underlying continuum from a fit to the 
/im region and not as we did from fitting a straight line between 8 and 13 /im. 



^In contrast to our analysis I Green et all (|2006i ) found the feature of V1057 Cyg to be comparable in shape 
(and thus in dust composition) to those of Bran 76, FU Ori, and V1515 Cyg. 
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object is the only one showing strong PAH emission bands. The three epochs of data for 
V1647 Ori allow us to study the variability of this object in the 10 fim region over a period 
of approximately five months. Figure [TT] shows that between October 2004 (epoch 1) and 
mid of March 2005 (epoch 2) the flux level decreased significantly. At the third epoch (end 
of March 2005) the flux appears to have slightly increased again (see also Figure [5]). The 
overall shape of the feature during the six months period did, however, not change. 

To put these results in a broader context with other young objects and to get a first 
idea on the dust grain properties, we plot in Figure [T^ the flux ratio at 11.3 and 9.8 /im 
against the computed peak over the continuum in the normalized spectra, i.e., the maxima 
of equation (1). This figure also shows the region typically occu pied by young TTauri 



stars and the slightly more massive Herbig Ae/Be sta. rs (see, e.g., iPrzygodda et al.l 12003 
van Boekel et al.ll2003l . l2005l : iKessler-Silacci et al.ll2006l ). We find that the objects presented 
here tend to have in general a relatively weak peak over continuum emission and a flux ratio 
between 0.8 and 1.0. While the strength of the peak over the continuum is interpreted as 
a tracer for grain sizes (with higher peak values denoting smaller grains), the flux ratio is a 
more general tracer for grain processing, i.e., grain growth as well as crystallization. This 
is explained by the fact that not only t he growing of grains le ads to a broader and flatter 
silicate feature between 9.8 and 12.0/im (IBouwman et al.ll200ll ). but also the crystallization 
proce ss introduces distinct em ission peaks of forsterite and enstatite longward of 10 fim (see. 



e.g., Ivan Boekel et al.ll2005l ). Taking these considerations into account, we find that our 



sample shows clear evidence for grain growth and dust processing. For a more quantitative 
analysis of the dust composition we fitted a dust model to the data as explained in the 
following subsection. 



3.2.2. Dust composition 

In Figure [13] we fitted an analytical dust model to the spectrum of FU Ori between 
7 and 17/im. To equally weight the low- and high-resolution part of the spectrum, we re- 
binned the data to a spectral resolution of 100, still sufficient to disentangle contributions 
from different dust species. Instead of fitting a two component model to the data, with one 
temperat ure for the underlyi ng continuum and one temperature for the emission layer, as 



we did in lOuanz et al.l (120061 ). we use a new model described in Juhasz et al. 2007 (to be 
submitted). In this model the temperatures for the disk midplane, the disk surface layer and 
an inner disk rim are fitted by powerlaw distributions instead of single black body temper- 
atures. In this respect, the model represents more realistically the temperature distribution 
in a protoplanetary disk. The input parameters for the model were taken from the simple 
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accretion disk model in iQuanz et al.l (120061 ) . Table H] summarizes the fitted parameters of 
the individual temperature components. 

The dust model was based on opacities for six dust species (amorphous olivine and 
pyroxene, crystalline forsterite and enstatite, silica, and amorphous carbon) in three grain 
sizes (0.1, 1.5, and 6.0/im). References for the optical p r operti es of the grains are given 
in Juhasz et al. 2007 (to be su bmitted) and iQuanz et al.l (120061 1 . In contrast the analysis 
presented in buanz et al.]j2006h . this time we included also carbon grains in the fit. However, 
as shown in Figure [13] and summarized in Table [5], the best fit (x^ = 45.3) does not require 
any carbon grains. The derived mass fraction s agree well with t he results based on the 
more simple two component dust model used in lQuanz et al.l (120061 ). While a large fraction 
of big, amorphous grains contribute to a high-temperature continuum, the smaller grains 
are responsible for the observed emission feature. Still, we emphasize that, according to 
our analysis, also in the disk surface layer larg e grains are present: Comparing the derive d 
composition to that of the interstellar medium (IKemper et al.ll2004l : lvan Boekel et al.ll2005l ). 
we find a significantly higher fraction of 1.5 /im grains in the dust of FU Ori. Furthermore, 
the model does predict a negligible amount of crystalline grains (< 0.01), no silica dust and, 
as mentioned above, no carbon grains. Certainly, the values for the derived mass fractions 
should not be taken literally, as, for instance, the temperature structure in the disk and 
the dust opacities are not perfectly known. However, the models presented here and in 
Quanz et al.l (120061 ) agree on two very important aspects: (1) there is clear evidence for 
dust coagulation in the spectra of FUors, and (2) there is hardly any contribution from 
crystalline grain s detectable in the d ata. This lack of crystalline grains was already noticed 
and discussed in iQuanz et al.l (120061 ) . but we will reinvestigate this issue below. 



3.2.3. The PAH features of Parsamian 21 



As its features are significantly different than those of the other FUors, we plot in 
Figure [H] the dereddened spectrum of Parsamian 21 and analyze it in more detail. The 
vertical lines illustrate, where typically PAH emission bands are located. While prominent 
PAH emission is present at 6.3, 8.2 and (possibly) 11.3 /im, no bands are seen at 7.7, 8.6 
and 12.7/im. In addition, the spectrum bears signs of underlying silicate emission and 
it is difficult to disentangle in the 11.3 /xm emission band the possible contributions from 
crystalline forsterite and PAHs. 



PAHs are commoii ly seen in the spectra of HAeBe stars (e.g.. lMeeus et al.ll200ll : lAcke &: van den Ancke 



2004 



Sloan et a. 



tra (iGeers et al 



20051) and, more recently, PAH emission was also detected in TTauri spec- 
2006h . However, the spectrum of Parsamian 21 differs significantly from the 
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PAH spectra of other young objects: it shows a broad and strong fe ature at 8.2 /xm whic h 
is barely detected and the main characteristic of group C objects in iPeeters et aLl ( 120021 ) . 
These objects show no PAH complex at 7.7 /x m and S.GjUni but the peculiar broad band 
at 8.2 /im whose origin yet remains unknown (IPeeters et al.ll2002l ). In section we will 
critically review the findings presented here and discuss their implications with respect to 
the FUor status of Parsamian 21. 



3.2.4- Qualitative analyses of the d-Sjim region 



As already noted and explained by lGreen et al.l (120061 ) FU Ori, Bran 76, V1515 Cyg, and 



also V1057 Cyg show absorption bands of gaseous H2O at 5.8 and 6.8 /im (Figure [5]) coming 
from a collection of rotation-vibration bands|l| V1647 Ori, however, shows evidence for ice 
absorption bands at 6.0 and 6.85 /im (see also Figure [20]) . We mention already here that 
this objects shows an additional ice feature at 15.2/im due to CO2 (Figure [T^ . V1647 Ori 
is the only object in the sample showing the 10 /im silicate feature in emission accompanied 
by significant ice absorption bands. A more detailed analysis of these ice features is given 
in section [3l3l Finally, as we will discuss in the following section, Parsamian 21 shows clear 
evidence for a PAH emission feature at 6.2 /im. 



3.3. Objects with 10 jim absorption 

3.3.1. Analyses of the lOfim region 

In Figure [13 we plot the optical depth in the 10 /xm region for the objects showing the 
silicate band in absorption. To derive the optical depth, we fitted the continuum with a 
polynomial of first order to the observed spectra between 8 and 13 /im, assuming that no 
absorption due to silicates is present at these wavelengths. Also here the exact shape of the 
continuum is difficult to estimate and any fit always carries uncertainties. However, even fits 
with higher order polynomials did not alter the results significantly and the main conclusions 
remained valid. Based on the assumed continuum we then computed the optical depth r 
using 

Fobs(A) = Fcont(A) -e-" (2) 



^Interestingly, already Sato et al. ( 1992f) found evidence for water vapor absorption between 1 and 3/im 
in the spectra of FU Ori, V1515 Cyg, and V1057 Cyg, but also L1551 IRS 5, V1057 Cyg, and Z CMa. 
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where -Fobs(A) is the observed flux and -Fcont(A) the assumed continuum. To derive the 
wavelength position and value of the maximum optical depth, we did the following: As some 
spectra are quite noisy (e.g., V1735 Cyg, RNO IB, Z CMa, and Reipurth 50) we fltted 
the data between 8.8 and 10.2 /im, i.e., the region with the greatest optical depth, with a 
polynomial of second order (shown as red, dash-dotted lines in Figure [15]). The wavelength 
position of the maximum optical depth (vertical, red, dash-dotted hues in Figure [15]) and 
the corresponding value are then derived from the flt. From the spectra with a high signal 
to noise ratio and strong absorption proflles (e.g., the Spitzer spectra of L1551 IRS 5 and 
HL Tau) it becomes clear that the flt matches the shape of the absorption feature quite well 
and that this approach seems reasonable. The black, dotted horizontal line in Figure [15] 
shows the assumed continuum level. In Table [6] we summarize the derived maximum optical 
depth for each object TsiUcate and the corresponding wavelength. We give also values for the 
extinction in the silicate band Agjiicate and estimates for the corresponding extinction in the 
optical Ay. Asiiicate was derived via 



A. 



silicate 



-2.51ogio (e" 



''"silicate 



(3) 



The errors for TsiHcate and AgiUcate are based on the 1-a uncertainties in the spectra and the 



resulting fluctuations of the polynomia l flt. Ay and its error are de r ived 



the extinction models of iMathisI (ll990l ) and IWeingartner fc Drainel ( 12001 
differe nt values of Ry for each model (Ry = 3.1 and 5.0 for lMathisI (ll990[) and Ry 



rom averagmg over 
) and assuming two 
3.0 and 



5.5 for IWeingartner fc Draind (l200l[ )). As already noted by I Cohen fc Kuhil (|l979l ) one has 
to be careful when deriving a value for the optical extinction based on the silicate feature, 
as uncertainties in the underlying continuum and unknown dust compositions influence the 
results. Hence, for comparison. Table [6] gives also values for Ay found in the literature. 
One has to keep in mind, though, that also here certain assumptions have been made and 
partly different observing techniques were apphed potentially leading to discrepancies in the 
derived values for Ay. 

For V346 Nor and RNO IC the values for Ay agree quite well. Also for Z CMa, Reipurth 
50 and L1551 IRS 5 the values are in general agreement. Due to the high extinction toward 
the latter two objects, a good estimate of Ay is difficult. W e would like to po i nt ou t, however. 



that a value of Ay > 15 mag for L1551 IRS 5 as found by Campbell et al.l (Il988l) and men- 
tioned, for instance, bv Istocke et~aD Jl988h . Iwhite et al.l J200oh and iFridlund et aP J2002h 
appears to be an overestimate. In yet unpublished NIR data taken with the Omega2000 
camera at the Calar Alto Observatory (Spain) we do see a near-infrared (NIR) counterpart 
to the L1551 IRS 5 system with Kg?5i9.7 mag. Even if this was purely scattered light from 
a disk the intrinsic K magnitude of the central object (s) (it might be a multiple system) can 
not be lower. Thus, the observed K-magnitude is a lower limit for the intrinsic K-luminosity 
of the source. The above mentioned value for Ay would translate into Ak ~ 15 mag, and 
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thus the absolute K-band magnitude of L1551 IRS 5 would be ~ — 11 mag for an assumed 
distance of 140 pc. For young low-mass objects this appears to be definitely too luminous 
even if accretion luminosities are taken into account as well. For V1735 Cyg and HL Tau 
our derived values for Ay are lower than those found in the literature. However, for VI 735 
Cyg no errors are given for the value found in the li terature. For HL T a u the optical ex- 
tinction derived from the radiative transfer model by iMen'shchikov et al.l (119991 ) may suffer 
from insufficient information about the source geometry or dust opacities which can easily 
lead to large uncertainties in the computed figures. The comparison for these objects clearly 
shows that, in particular for embedded objects, it is difficult to derive consistent values for 
Ay if different techniques are applied. 

Finally, RNO IB a lso shows a discrepancy between our value for Ay and that given by 
Staude &: Neckell (Il99ll ). In this case, however, the silicate absorption feature bears evidence 
for additional superposed emission as will be described in the next paragraph. Thus, we 
attribute at least part of the missing optical depth to an underlying silicate emission feature. 

To measure not only the depth of the absorption profile, but to get also an idea of 
the dust composition responsible for the observed absorption, we analyzed the shape of the 
silicate feature. In Figure [16] we show the optical depths computed or measured for different 
dust grain compositions. Depending on the composition the wavelengths of the maximum 
optical depth changes. The most extreme cases are amorphous olivine grains with a size of 
0.1 nm peaking at a wavelength slightly longward of 9.7/im, and amorphous pyroxene grains 
with the same size peaking around 9.2 /xm wavelength. As mentioned above, in Table [H] we 
already give the wavelengths of the maximum optical depths observed toward our sources and 
it shows that there is indeed a broad range of values, indicating different dust compositions. 

In Figure [T7] we plot again the observed optical depth of our objects as shown in Fig- 
ure but now we overplot one or two of the absorption profiles presented in Figure [121 
scaled to the maximum observed optical depth. Whether this scaling is appropriate is diffi- 
cult to determine but it shows that most observed features agree reasonably well with one or 
two of the reference profiles. Three objects (L1551 IRS 5, V346 Nor and V1735 Cyg) show 
additional absorption longward of llyum when compared to the dust fe atures, but at least 



part of this absorption can be attributed to H2O libration bands (e.g., iPontoppidan et al. 
2005! ) . VI 735 Cyg, RN O IC and Z CMa seem to agree best with the dust model com- 
puted by iDraind (120031 ) based on the particle size distribution from IWeingartner fc Draine 
(120011 ). V346 Nor requires shortward of 10 /xm a mixture of the sarn e model and the dust 
composition observed toward the galactic center (IKemper et al.ll2004l ). Longward of 10 /im 
the former model alone fits better. For Reipurth 50 a combination of bo th models explains 



the data well over the whole wavelength regime with the iDraind (120031 ) model fitting the 
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short and the iKemper et al.l (l2004j ) model fitting the long wavelength part. L1551 I RS 5 is 
best fi tted with a mixture of small, amorphous pyroxene grains and the model from iDraine 
(12003). The profile o f HL Tau agrees extremely well with the galactic center dust profile from 
Kemper et al.l (120041 ). To account for a small shift toward shorter wavelengths a little bit 



more amorphous pyroxene seems to be required or small uncertainties in the position of the 
peak absorption from the continuum fit can explain this shift. Finally, the absorption profile 
of RNO IB is special in two ways: First, the short wavelength range is best fitted solely 
with small, amorphous pyroxene grains. And secondly, between 10 and 11.5 /im the profile 
shows a significantly lower optical depth than any of the reference profiles. In combination 
with our findings for the visual extinction (see above) this suggests that RNO IB shows 
silicate emission superposed on the absorption feature. Interestingly, the apparent decrease 
in optical depth longward of 10 /xm indicates that at these wavelengths the contribution of 
the emission is relatively larger compared to the other wavelengths in the silicate band. This 
in turn means that the emission profile does not have the typical shape of the ISM dust 
feature but a broader and more evolved profile like the other FUors emission profiles shown 
in Figure [TOl 



3.3.2. The 15.2 CO2 ice feature 



The shape of the bending mode profile of CO2 ice around 15.2 /im does not only allow a 
detailed analysis of the involved ice inventory, but it enables us to derive in formation about 
potential ice processing due to heating effects (e.g.. |Pontoppidan et al.ll2005l ). In Figure [T8] we 
present two reference cases which we will use in the following as a benchmark for comparison 
with the FUor spectra. The first obj ect in Figure [TSl (0 X2) is a highly extincted background 
star behind the Serpens dark c loud (IKnez et al.1 120051) . and the second object (IIII46) is an 
embedded low- mass protostar (iBoogert et al.ll2004l ). 



While the spectrum of CK2 probes mainly the ice composition of the intervening dark 
cloud, the spectrum of HH46 bears information about the immediate environment of the 
young star. Apart from being slightly broader, the spectrum of IIII46 shows a double- 
peaked sub-structure which is caused by crystallization and effective segregation of the CO2 
and H2O ice involved in the absorption. Thes e processes take pla ce when ice mixtures 
with concentrations of C02/H20>1 are heated (IBoogert et al.ll2004l ). In space, this phase 
transition from amorphous to crystalline ice is expected to occur between 50-90 K and 
thus at higher temperatures than typically found in cold, dark molecular clouds. Thus, 
the 15.2 /xm ice feature of HII46 shows that (at least) part of t he ic e must already have 
been heated by the embedded protostar. In fact, IBoogert et al.l ( 120041 ) fit the feature with 
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a two-component ice model based on laboratory spectra with one component being highly 
processed polar ice with a laboratory temperature of Ti^b = 125 K, and the other component 
being an H^O-rich, CH.sOH -deficient cold ice with Tiab = 10 k|^ The spectrum of CK2 was 



fitted by iKnez et al.l (120051 ) solely with cold ice components. They used a polar mixture 
of H20:C02=1:1 and H2O:CO2=10:l at 10 K with a ratio of 2:1 and an additional apolar 
component of CO:N2:CO2=100:50:20 at 30 K. The overall polar fraction was assumed to be 
78%. 

In Figure [19] we present the 15.2 /xm features observed toward our FUor sample. To 
derive the optical depths we fitted the continuum with a straight line fixed around 14.65 
and 16.3 /xm. For each object we overplot either the (scaled) spectrum of HH46 or CK2 
depending on the shape of the profile. While the profiles of V1647 Ori and V346 Nor agree 
better with the profile of CK2 (the spectrum representing unprocessed ice), L1551 IRS 5, 
RNO IB and IC, and HL Tau show evidence for a double-peaked sub-structure and thus 
heating effects and processed ice. The comparatively bad quality of the VI 735 Cyg spectrum 
does not allow a solid comparison to either reference spectrum. It is noteworthy that the 
spectrum of V1647 Ori is almost an exact copy of CK2 indicating that the ice composition is 
mostly identical. Based on the fitted reference spectrum we computed also the optical depth 
Ti5.2Atm for each object and summarized the results in Table [71 The errors are based on the 
l-cr uncertainties in the observed spectra. 



3.3.3. The 6.0 and 6.8 nm ice features 



Although frequently observed toward high- and low-mass sources, the two well-known 
ice features at 6.0 and 6.8 /zm are quite complex and difficult to interpret. Certainly, a large 
fraction of the optical depth of the 6.0 /im band ca n be attributed to H2O ice, but also other 
species might contribute to this absorption feature (iKeane et al.ll200ll ). For instance, slightly 
short ward, at roughly 5.85 Atm, an additional ab sorption shoulder is sometimes superposed 



(e.g.. lPontoppidan et al.ll2005l : iKeane et al.ll200ll ) for which formaldehyde (H2CO) and formic 
acid (HCOOH) are theoretical candidates. Without any additional information (e.g., the 
3.08 /im band of H2O or the 3.47 /xm band of H2CO) it is thus difficult to determine the true 
water ice content in the 6.0 /xm band. 

An additiona l absorption feature at 6.85 /im is often observed towar d protostars (e.g. , 
Keane et al.ll200ll ). but also toward the extincted background star CK2 (IKnez et al.ll2005l ) 



^The presence of CH3OH in the ice feature can potentially be traced by the shape of the long-wavelength 
wing of the CO2 profile showing additional absorption if CH3OH is present in higher abundances. 
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and the edge-on disk CRBR 2422.8-3423 flPontoppidan et al.lboosh . Although a final identi- 
fication of this band has yet to be provided, NHj" seems to b e one of the most promising can- 
didates (ISchutte fc Khannall2003l : iPontoppidan et al.ll2005l ). However, Ivan Dishoeckl (120041 ) 
for example, mention also methanol (CH3OH) as a potential carrier of this absorption band. 

Given all these uncertainties we restrict ourselves in this section to the computation of 
the optical depth of both of the above mentioned absorption bands. Like in the previous 
section we assumed a straight line for the continuum anchored at 5.4 and 7.6 /zm. To derive 
the optical depths we then fitted a polynomial of fourth order to both absorption dips to 
eliminate the noise in the spectra. Figure [20] shows the observed spectra between 5.5 and 
7.5 fim on an optical depth scale (black lines) with the resulting fits overp lotted (red, dashed- 
dotted lines). The computed optical depth for each band and each object is given in Table [71 
Like in the previous section the errors are based on the l-cx uncertainties in the observed 
spectra. 



3.4. The fading of OO Ser and V1647 Ori 

For two objects (00 Ser and V1647 Ori) we have multi-epoch data and can derive 
some conclusions on the variability of these objects. As illustrated in Figure [U we fitted 
a straight line to the data of 00 Ser between 15 and 30 fim to estimate the decay in fiux 
density observed over the 5 epochs (see, Tabl^I])- The wavelength range was chosen as in 
this regime all spectra are still relatively clear of artefacts and spikes partly seen at longer 
wavelengths. In Table M we summarize the fiux density level at 20 and 30 fim at each epoch 
and give also the slope measured between 15 and 30 /im. The errors are derived from the 
goodness of the fit to the data. Between the first observations on April 14, 1996, and the last 
observations, carried out September 22, 1997, the fiux densities decreased to roughly 50% of 
the initial values. These data demonstrate that 00 Ser faded rapidly over relatively short 
timescales and that it might be an intermediate object between a typical FUor and an EXor 
having fading timescales of several d ecades and months, respectively. Based on photometric 



monitoring at infrared wavelengths, iKospal et al.l (l2007l ) came to a similar conclusion and 



predicted that 00 Ser should return to its pre-outburst luminosity not before 2011. 

In Table [8] we also summarize the fiux densities of V1647 Ori observed at 8, 20 and 
30 fim at three different epochs. The errors are taken directly from the Spitzer spectra. It 
shows that also this object faded significantly over a period of less than 5 months between 
October 2004 and March 2005. In addition, short time variations in the fiux levels seem to be 
present, as between the two epochs in March 2005 the object became slightly brighter again. 
In consequence, these data support the assumption that the outburst of V1647 Ori may also 



be intermediate between FUor- and EXor-type events similar to 00 Ser (iMuzeroUe et al. 



20051 : lAcosta-Pulido et al.ll2007f ). 



3.5. Additional emission lines 



For completeness we show in Figure [21] absorption and emission lines identified in the 
high-resolution regim e of the Spitzer sp ectra, part of which are difficult to identify in Fig- 
ures [7] and [HI Already [Green et al.l (|2006[ ) noted the [S III] emission lines at 18.7 and 33.4 /xm 
in the spectrum of V1515 Cyg and argued that they originate from extended emission in the 
region and not from the object itself. The spectra of RNO IB and RNO IC show evidence 
for H2 quadrupole emission around 17/im a nd additional lin es in the low-resolution part 
of the spectrum at shorter wavelengths. In [Quanz et al.l (|2007[ ) we analyzed these emission 
lines in detail and concluded that they are related to shocks within a molecular outflow 
powered by the nearby embedded object IRAS 00338-1-6312. [Fe II] lines around 17.9 /xm are 
present in the spectra of L1551 IRS 5 and RNO IB, and L1551 IRS 5 shows also the [Fe II] 
line near 26.0 /xm. The [Fe II] lines of L1551 IRS 5 were already detected in the ISO/SWS 
spectrum of this source and att ributed to hot and dense material located close to the root of 
the outflow ( IWhite et al.[[2000[ ). The line intensities did not fit the predictions from shock 
models. Concerning the [Fe II] line in the spectrum of RNO IB, it see ms likely that i t arise s 
in the outflow shocks that also excite the H2 emission lines discussed in [Quanz et al.[ (l2007[ ). 
Finally, V1057 Cyg shows gaseous CO2 absorption slightly shortward of 15.0 fim. 



4. Discussion 



4.1. Two categories of FUors 



The results presented here suggest that the sample of FUors can be divided into two 
categories based on the observational appearance of the 10 /im silicate feature. We decided 
to call the absorption feature objects "Category l"-FUors (9 objects in Table [3]) and the 
emission feature objects "Category 2"-FUors (6 objects in Table [3]). In the following we 
discuss the properties of the two categories in more detail. 
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4-1.1. Category 1 FUors: silicate and ice absorption features 
In general it is possible to observe the silicate feature in absorption if the circurnstellar 



disk surrounding a young object is seen close to edge-on. iMen'shchikov &: Hennind (119971 ) 
showed that typically for a disk with an opening angle of 20° between the upper and lower 
disk surface the silicate feature appears in absorption only if the disk inclination is < 10° from 
edge-on. Given the amount of objects showing silicate absorption compared to the number 
of objects with silicate emission and assuming a random distribution of the orientation of 
the accretion disks, it is unlikely that all Category 1 FUors are seen edge-on. Rather, 
these objects are still more deeply embedded in their molecular envelopes covering a larger 
solid angle than the edge-on disk alone. These envelopes cause the ice and dust absorption 
features. 

Figure [17] illustrates that all but one silicate absorption bands are fitted best either 



with the silicate composition observed t oward the Galactic cen t er by iKemper et al. (120 04*) 
or with the astronomical silicates from IWeingartner &: Draind ( 200ll ) and iDraind (2003|)|^' 



This finding shows that all absorption features can be fitted with small, amorphous silicates 
and hence the extinction is caused by pristine and not processed dust. Only the spectrum of 
RNO IB is difficult to fit with any ISM dust compostition, but, as already mentioned, there 
are strong hints that the absorption feature is altered by a superposed emission feature. 

In Figure [22] we show the observed optical depths of the ices and the silicate feature in 
a scatter plot to search for any correlation. Although the range of optical depths we probe 
here is limited, it seems that at least for the ices there seem to be correlations (left and right 
plot in Figure [22]). Fitting a straight line to the data we find: 

ri5.2^m = (0.349 ± 0.037) • T^,o,.rn + (0.040 ± 0.009) (4) 

T-6.8,.m = (0.858 ± 0.067) ■ rg.o^n, + (0.007 ± 0.016) (5) 

This implies that the physical and chemical conditions within the envelopes and clouds 
causing the absorption are similar. A more detailed analysis, e.g., the determination of 
different ice abundances relative to water ice, is beyond the scope of this work. For such 
an analysis the water ice feature around 3 /im is required as it suffers less from additional 
contributions of other ice species compare to the 6 /im feature. 

For a correlation between the optical depths of the silicates and ices the situation is 
slightly different because, as we have already pointed out, the silicate feature can be influ- 



'L1551 IRS 5 requires a fraction of additional amorphous pyroxene grains. 
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enced by superposed emission. This shows nicely in the middle plot in Figure [22|, where 
RNO IB is shifted with respect to the other objects. 

Finally, it is interesting to note that the objects showing evidence for ice processing 
in Figure [19] tend to show higher optical depths in the ice features than the other sources. 
This might suggest that the extinction for the latter objects (V346 Nor, V1735 Cyg) might 
be caused by ices somewhere in the line of sight to the source, rather than by material 
related to the young star. This might also explain why V1647 Ori shows weak silicate 
emission, but ice absorption: the extinction is caused by cold foreground material, reflected 
also in the observed high value for Ay (see caption Figure [TOj). Ground-ba sed observations 



find the spectrum to be fiat in the 10 /xm regime (lAbraham et al.l l2006l ). However, the 
sensitivity of Spitzer allows the detection of a weak feature resulting from an extinguished 
but intrinsically strong sihcate emission band. 



4-1-2. Category 2 FUors: the silicate emission feature and its dust composition 



As presented in Figure [TO] and described in the related section, the emission profiles 
bear evidence for dust grain processing. Even after the correction for apparent interstellar 
extinction, the shape of the silicate profile differs from that of typical ISM dust. Like for 
TTauri stars and the slightly more massive HAeBes, the origin of the emission feature in 
the spectra of the FUors is the heated surface layer of the accretion disk. The apparent 
grain processing is believed to be only possible in circumstellar disks and not in he less dense 
circumstellar envelopes. While the emission layer for TTauri stars and HAeBes is mainly 
heated be the central stellar object, for FUors the hot inner par ts of the accretion disk itself 
can act as the main illuminating source (e.g., iLachaumd l2004j ) . Due to the high accretion 
rates of FUors those inner regions are extrem ely hot and accoun t for a significant fraction 
of the total fiux even at optical wavelengths (iQuanz et al.ll2006l ). Furthermore, accretion 



disk mod els with a flared g eometry were not only able to explain the SEDs of FU Ori and 
Bran 76 ( Green et al. 20061) . but they co uld also reproduce interferometric observations in 



the NIR (iMalbet et al.l l2005h and MIR (lOuanz et al.1120061 ). V1057 Cyg and V1515 Cyg 
show more emission at longer infrared wavelengths than the previous objects, which can be 
accounted for assuming a remnant infalling envelope in a ddition to the accretion dis ks. This 
is supported by the results from K-band interferometry iMillan-Gabet et al.l (120061 ) , where 
envelopes are required to explain the low NIR visibilities. 

Concerning the dust composition it is interesting to note that, given the low peak- 
over-continuum ratio illustrated in Figure [121 grain growth must have already set in. This 
is supported by the results of our dust model fit. Furthermore, the spectra do not show 
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evidence for crystalline dust particles. The mass fraction of crystalline particles in the dust 
model computed for the spectrum of FU Ori was negligible (see, Table [5]). Since strength 
and shape of the silicate feature of FU Ori is comparable to those of Bran 76 and V1515 
Cyg (Figur jTTl) . the dust composition in all objects is similar. The spectrum of V1057 Ori 
is even broader and less pronounced than the other spectra, indicating even larger grains in 
the disk surface layer. The intrinsic feature of V1647 Ori is stronger than those of the other 
FUors, but no prominent signs of crystalline silicates are present either 



As already mentioned in lQuanz et al.l (l2006l ). there are several reasons for which stronger 
crystalline features could have been expected to be detected. The high disk accretion rates 
should ease the detection o f crystalline part icles in two ways: (1) high accretion rates lead 
to high disk temperatures (IBell et al.l 119971 ) which in turn should increase the amount of 
crystalline particles produced by annealing processes at T>800 K. (2) an incr ease in the 
accretion rate should also increase the radial and vertical mixing in the disk (IGaill l200ll ) 
transporting the crystalline particles farther out and to the disk surface, where they can 
be detected by means of MIR spe ctroscopy. Howeve r, neither in the innermost disk regions 
probed with MIR interferometry (iQuanz et al.ll2006l ) nor in the Spitzer spectra presented 
here in Figure (TU] or at longer wavelengths in Figure [7] we see any striking evidence for 
crystallinity. This means, that either those grains do not exist in large amount in these disks 
or they are somehow hidden. 

The disks of FUor objects are different from those of TTauri or Herbig star disks in 
the sense that in the radially innermost regions of FUor disks the luminosity is accretion 
dominated (TTauri and Herbig star disks are irradiation dominated everywhere). At larger 
radii, also FUor disks are irradiation dominated, albeit the main heating source may not be 
the central star itself, but rather the hot inner disk regions close to the star where accretion 
dominates (e.g.. lLachaumell2004l ) . Disk regions that are irradiation dominated have a surface 
layer that is warmer than the underlying disk interior; such an " inverted" temperature profile 
causes dust features to appear in emission. In the accretion dominated regions, the main 
heating source is the release of gravitational energy in the disk midplane. Here, the disk 
interior is at least as hot (and likely hotter) than the disk surface, causing dust features to 
be effectively hidden, or even to appear in absorptior0. Thus, it is conceivable that there 
are significant amounts of crystalline silicates present in the hot inner disk regions of our 
targets, even if they do not show up prominently in the spectra. 



^"^The spectrum of Parsamian 21 does show a prominent feature at 11.3/im but we attribute most of the 
related flux to PAH emission and not to crystalline silicates (see section [3. 2.3p . 

^^The gas absorption features between 5 and 8 /zm described in 3.2.4 likely originate in this hot disk region. 
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But what about crystalline silicates further from the central object, in irradiation domi- 
nated regions? Given the high accretion luminosity in our objects it seems unavoidable that 
large amounts of crystalline silicates are produced in the innermost disk reg ions. In addition 



to th is, crystalline silicates may be produced in situ at large radii in shocks (IHarker &: Desch 



2OO2I ) or electric discharges (jPilipp et al.lll998l : iDesch fc Cuzzill200Cll ) . If the disk is well mixed 



in the vertical direction, the crystalline silicates should be present in the surface layer of the 
disk at large (^few AU) radii, and show up in emission, which they evidently do not (at 
least, to a much lesser extent than in many TTauri and Herbig Ae/Be star disks, which 
are more evolved). This inevitably leads to the conclusion that crystalline silicates are not 
abundant in the surface layer of FUor disks at radii of more than a few AU from the central 
object, contrary to expectations. 

At this point, we have no proper explanation for this observation. In order to animate 
the discussion, we will post the following idea, but stress that at this point this is mere 
speculation. The crystalline silicates are not present in the surface layer at larger radii, but 
may be present in the disk interior, where they are expected to be formed by virtue of the 
processes mentioned before. Thus, we need to explain why the crystalline silicates do not 
get mixed in the vertical direction all the way to the disk surface. One possibility is that 
the crystalline particles are somehow less well coupled to the gas, possibly because they are 
more compact (or less "fluffy") than the unprocessed dust particles. This may cause the 
selective settling of crystalline dust. A second possibility is the following. The accretion rate 
in FUor objects is high, suggesting that the disk is still being supplied with fresh material 
from the maternal envelope. This material, which is thought to contain exclusively small, 
amorphous dust particles, will in part "rain" onto the disk surface at larger radii. If the 
rate at which this material falls onto the disk is higher than the rate at which it is mixed 
through the disk in the vertical direction, the disk may be covered by a "blanket" of pristine 
material, effectively hiding the more processed material present in the disk interior. It is 
unclear whether this idea can be harmonized with the observation that the dust we see at 
the disk surface has undergone significant grain growth. This issue clearly calls for a much 
more detailed and quantitative investigation, which is beyond the scope of the present work. 



4.2. Unifying the two categories of FUors: An evolutionary sequence 

The spectral properties of the sample presented here draw a rather inhomogeneous 
picture of the group of FU Orionis objects. However, as most of them are convincingly 
classified as FUors, they do share some common observational and physical properties. In 
the following, we present an idea, how the apparent differences can be explained within a 
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unified paradigm of FU Orionis objects. 



Already 20 years ago iHerbigl (119771 ) and iHartmann &: KenyonI (Il985l ) suggested that 
each FUor might undergo several outbursts. More recent theoretical models show that 
gravitational instabilities in the accretion disk, driven by continuing infall from a remnant 
envelope, can account for intense bursts of high accretion rates which a r e intersecting mor e 
prolonged, quiescent periods of low accretion (IVorobyov fc Basul l2006l : iBoley et al.l l2006l ). 
After several of these outbursts have occurred within a time span of several 10^ years, the 
envelope, which is the trigger of the disk instabilities, vanishes, and the object enters finally 
a state of permanently low accretion. If one relates this final phase to the classical TTauri 
phase of low-r nass YSOs, then FUors are younger than classica l TTau ri stars as already 
pointed out by IWeintraub et al.l (Il99ll ) and lSandell &: Weintraubl (120011 ). 



Kessler-Silacci et al.l ( 120051 ) analyzed the 10 /xm silicate feature of 34 young stars and 



confirmed that the evolution of the feature of low-mass stars and the overall SEDs is similar 
to that of intermediate-mass stars: embedded objects, showing a pristine silicate band in 
absorption, evolve into objects showing a combination of silicate absorption and emission, 
and finally pure emission features appear, where dust processing leads to a broad range of 
shapes and strength. Since all these stages are represented in our sample, we believe that, 
indeed, the FUor-phase is the link between the embedded Class I objects and the more 
evolved Class II objects. In this context, the objects showing silicate absorption are younger 
and at the beginning of the period where subsequently numerous FUor outbursts will occur, 
while the objects showing pure emission features are more evolved and possibly near the end 
of their FUor period. As illustrated in Figure [23] the s pectral indices of t he emission objects 
are indistinguishable from Class II objects in Taurus (IFurlan et al.l 120061 ). In particular, the 
objects FU Ori and Bran 76 appear to be the most evolved objects, because, as mentioned 
above, their SEDs do not require the presence of a large remnant envelope. A nice example 
of an intermediate object is RNO IB, where the dominating absorption feature is altered by 
the underlying silicate emission from the accretion disk. Figure [2l] illustrates with simple 
sketches the main features of the the two categories of FUors. It should be mentioned that 
theoretically a Class II FUor seen close to edge-on might be interpreted as a Class I FUor. 
As outlined in 4.1.1 it is, however, statistically very unlikely that all of the Class I objects 
shown here are Class II FUors "in disguise". 

One question arises, though. If FUors are younger than classical TTauri stars but do 
already show clear evidence for dust coagulation, then why do some classical TTauri stars 
show silicate emission features consisting purely of pristine and unprocessed dust? One 
possible explanation is linked to the high accretion rates of the FUors. It is believed that 
high disk accretion rates are related to a higher degree of turbulence and thus mixing in the 
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disk. As a result, larger particles are coupled to the turbulent gas and mixed throughout the 
disk and to the disk surface, where they can be observed. Once the accretion rate, and hence 
the amount of mixing, drops, larger particles tend to settle to the disk midplane much faster 
than small particles, which remain in the disk surface layer. Thus, although large particles 
are present in the disk of a TTauri star, they might not reach the disk surface, and only 
the small grains produce a silicate emission feature. Again, as for the issue concerning the 
lack of crystalline grains, a more detailed theoretical modeling is certainly required to see to 
what extend this qualitative explanation is valid. 



4.3. Parsamian 21 - intermediate mass FUor or Post-AGB star? 



The strong PAH emission of Parsamian 21 is unique among the sample of FUors pre- 
sented here, and the analysis rise serious doubts w hether Parsamian 21 is indeed a member 
of this group. The initial classification was made by lStaude &: Neckell (119921 ) based on optical 
spectroscopy. These authors derived a spectral type of F5Iab, found a prominent P-Cygni 
profile in Ha and observed shock-induced emission in [O I], [N II] and [S II]. They also 
mention the detection of Li in their spectra, which, if present, does not seem to be signifi- 
cantly above the noise level in the data. New observations with higher spectral resolution 
and higher signal-to-noise would certainly be eligible. Although some of these features are 
commonly observed in FUors, they are not unique to this group of objects, as also HAeBes 
show, for instance, P-Cygni profiles in Ha. In addition, two of the main properties of 
FUors, an outburst in opti cal light and CO bandhead absorption profiles in the NIR (e.g., 
Hartmann fc Kenyorull996l ). have not yet been observed for Parsamian 21. Looking at the 
immediate surroundings of Parsamian 21 in optical data from the Digitized Sky Survey!^ it is 
also striking that no dark cloud complex is linked to Parsamian 21 and hence no connection 
to a larger star-forming region is present. 

To excite the PAH emission, sufficiently strong UV radiation is required. However, 
most FUors have a later spectral type t han Parsamian 21 and lack, unlike TTauri stars, in 
general UV continuum excess emission (IHartmann fc Kenyorull996l ). As mentioned above, 
the positions of the PAH emission bands are untypical for young st ars and so far onl y 7 
objects are known to show ch aracteristics of PAH class C objects (ISloan et al.l 120071 ) as 
defined by iPeeters et al.l (120021 ). Among these 7 objects only two presumably young stars 
(SU Aur and HD 135344) show the peculiar PAH band close t o 8.2 /xm. All other PAH class 
C objects are post-AGB stars or red giants (ISloan et al.l 120071 : iPeeters et al.ll2002l ). 



^^http: / /archive. eso.org/dss/dss 
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In addition, one of the Post-AGB stars from jPeeters et al.l (20021 ) and also HD 56126, 
a Post-AGB star with PAH and sihcate emission (IHonv et al.l 1200311. have the same spectra l 
type as Parsamian 21. However, we should mention that iDibail (119691 ) and lThe et al.l (119941 ) 
give a spectral type of A5Vea for Parsamian 21. This is too early for a typical FU Orionis 
object, but could also explain the PAH emission. Finally, the bipolar emission knots seen 
in Ha and [N II] in the optical spectrum of Parsamian 21 (jStaude fc Neckell Il992l ). are 
also found in the bip olar outflo ws of evolved stars and "pre-planetary" nebulae (e.g., in the 
"Butterfly Nebula", ISollI 120001 ). They are thus not only seen in young objects. 



Based on these findings, we believe that the FUor status of Parsamian 21 is at least very 
questionable. Either this object represents an intermediate mass FUor object, suggesting 
that also stars in this mass regime undergo phases of enhanced accretion, or, and this appears 
even more likely, Parsamian 21 is not even a young object but an evolved star, sharing typical 
properties with Post-AGB stars. 



5. Conclusions and future prospects 



Our conclusions can be summarized as follows: 



• We presented the first coherent space based spectroscopic MIR study of 14 FUors ob- 
served with the Spitzer Space Telescope or the Infrared Space Observatory. 
The sample includes roughly two thirds of the known FUors or FUor candidates. 

• Based on the appearance of the 10 /xm silicate feature, we divided the sample into two 
categories: Category 1 objects show the silicate feature in absorption, and the spectra 
show additional absorption bands at 6.0, 6.8 and 15.2 /im due to ices. Category 2 
objects show silicate emission and (most of them also) indications for water vapor 
absorption at shorter wavelengths. Only one Category 2 object (V1647 Ori) shows 
ice absorption bands which we explain by foreground extinction and an intrinsically 
strong 10 /im emission feature. 

• The silicate absorption is best explained with dusty and icy envelopes surrounding the 
Category 1 objects. Statistical reasons argue against all objects having accretion disks 
seen edge-on. The silicate emission of the Category 2 FUors arises from the surface 
layer of the surrounding accretion disks. 

• The shape of the silicate band of the Category 1 objects is in agreement with typical 
dust compositions of the ISM. For one object (RNO IB) the shape of the feature and 
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the decrease in optical depth longward of 10 //m can be explained with a superposed 
emission feature. 

• Optical depths for the observed sihcate and ice absorption bands were derived. We 
find an apparent correlation among the optical depths of the ices indicating similar 
environmental conditions for the objects. For the silicate feature, no correlation is 
expected. 

• Using different extinction curves we computed the optical extinction toward the objects 
based on the depth of the 10 //m silicate features. The results are in agreement with 
values found in the literature, given the uncertainties in the dust models and for the 
optical extinction values in the references. For RNO IB the derived value for Ay 
is smaller than expected, indicating that on-top silicate emission might influence the 
optical depth at 10 /im. 

• The emission profiles of the Category 2 objects show clear evidence for grain growth. 
Fitting a dust model to the spectrum of FU Ori reveals that, indeed, larger grains 
than typically observed in the ISM are required to explain the shape. As FUors are 
presumably younger than TTauri stars (see below), this indicates that grain growth 
sets in very early during disk evolution. 

• Despite the high accretion rates of the FUor accretion disks and the resulting higher 
disk temperatures and mixing rates, we find hardly any evidence for crystalline grains. 
So far we lack a clear explanation for this observational results and leave it to further 
investigations. 

• The two categories of FUors can be explained within a single paradigm, where Category 
1 objects are younger and similar to Class I objects, while Category 2 FUors are 
more evolved and show already properties of Class II sources. This explanation is in 
agreement with theoretical models which expect FUors to undergo several outbursts 
before they enter the more quiescent classical TTauri phase. Thus, the FUor-phase 
might indeed be the link between Class I and Class II objects and common to most 
young low-mass stars. 

• For 00 Ser and V1647 Ori the multi-epoch data allows an analysis of the post-outburst 
fading of the objects. Both objects fainted significantly over timescales of a few months, 
suggesting the outbursts might be intermediate between the long-lived FUor-eruptions 
and the short-lived EXor-type events. 

• Only one object (Parsamian 21) shows PAH emission similar to that often observed in 
Post-AGB stars. We find that most other observational data for Parsamian 21 can also 
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be explained with the object being an evolved star. In consequence, the FUor-status 
for this object is questioned. 

Based on these findings, future investigations might include the following points: 

• A complete MIR spectroscopic census of all known FUors would complement the data 
presented here and might help to derive conclusions concerning the duration of the 
FUor phase using statistical arguments. 

• Multidimensional radiative transfer models have not yet been applied to FUor accretion 
disk, but they are required to derive a coherent picture of the disk structure, including 
the emission layer. 

• Models for the dust evolution in accretion disks could try to explain quantitatively the 
observed large grain population and the apparent lack of crystalline silicates. 

• Finally, if our conclusions are correct, then some of the known low-mass Class I ob- 
jects might be FU Ori objects hidden in a quiescent phase between two consecutive 
outbursts. The observational properties of these objects might be revisited to search 
for any indication of FUor properties. 
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Table 1. Journal of ISO-SWS observations. The coordinates denote the pointing position of the telescope as saved 
in the header of the data files. The speed of the observations, the time spent on the target and possible pointing 
offsets in the telescope's y and z axis are given as well. The offsets were corrected for during the data reduction 

applying the measured beam profiles along the different axes. 



Object 


RA (J2000) 


DEC (J2000) 


AOT / Speed 


Time on 
target [sec] 


Offset ["] 
y / z 


Date 




00 Ser 


18h29m49.05s 


+01dl6'19.2" 


SWSOl / 1 


1062 


- /- 


1996-04- 


■14 




18h29m49.08s 


+01dl6'19.8" 


SWSOl / 1 


1140 


4/3 


1996-10- 


■24 




18h29m49.05s 


+01dl6'19.2" 


SWSOl / 2 


1140 


-/- 


1997-03- 


■08 




18h29m49.05s 


+01dl6'19.2" 


SWSOl / 3 


3454 


-/- 


1997-04- 


-12 




18h29m49.09s 


+01dl6'19.8" 


SWSOl / 3 


3454 


-/- 


1997-09- 


-22 


RNO IB 


00h36m46.24s 


+63d28'54.3" 


SWSOl / 2 


1912 


-/- 


1996-08- 


-27 


V346 Nor 


16h32m32.05s 


-44d55'28.9" 


SWSOl / 2 


1912 


-/- 


1996-08- 


■31 


Z CMa 


07h03m43.17s 


-lld33'06.6" 


SWSOl / 2 


3454 


3/3 


1997-11- 


-07 


V1735 Cyg 


21h47m20.60s 


+47d32'04.9" 


SWSOl / 2 


1912 


4/4 


1996-08- 


■06 


Reipurth 50 


05h40ml7.89s 


-07d27'29.3" 


SWSOl / 3 


3454 


5/5 


1997-10- 


■13 


L1551 IRS 5 


04h31m34.06s 


+18d08'04.8" 


SWSOl / 4 


6538 


-/- 


1997-09- 


■06 



Table 2. Journal of Spitzer-IRS observations. The coordinates denote the average sht position of the low- 
resolution spectrograph as computed by the on board software. The AOR of the observations, the time spent on the 
target for the different modules and the observation date are given as well. 



Obiect 


RA fJ2000') 


DEC fJ2000') 


AOR 


TntpEfration time fsecl 


Date 


Bran 76 (BBW 76) 


07h50m35.52s 


-33d06'24.12" 


3571200 


12 (SL, LL) 


2004-04-14 


FU Ori 


05h45m22.39s 


+09d04'12.5" 


3569920 


12 (SL, SH, LH) 


2004-03-04 


L1551 IRS 5 


04h31m34.08s 


+18d08'04.92" 


3531776 


6 (SL, SH, LH) 


2004-03-04 


Parsamian 21 (HBC 687) 


19h29m00.72s 


+09d38'47.11" 


5039872 


36 (SL) 

48 (SH, LH) 


2004-04-18 


RNO IB 


00h36m46.34s 


+63d28'53.76" 


6586624 


36 (SL) 

60 (SH, LH) 


2004-01-07 


RNO IC 


00h36m46.89s 


+63d28'58.44" 


6586624 


36 (SL) 

60 (SH, LH) 


2004-01-07 


V1057 Cyg 


20h58m53.76s 


+44dl5'28.44" 


3570176 


12 (SL, SH, LH) 


2003-12-15 


V1515 Cyg 


20h23m48.00s 


+42dl2'25.56" 


3570432 


12 (SL, SH, LH) 


2004-05-11 


V1647 Ori^ 


05h46ml3.13s 


-00d06'05.21" 


12261120 


12 (SL) 

24 (SH, LH) 


2004-10-20 




05h46ml3.15s 


-00d06'04.41" 


11569920 


48 (SL) 

484 (SH), 240 (LH) 


2005-03-11 




05h46ml3.14s 


-00d06'04.69" 


12644096 


12 (SL, LL) 


2005-03-24 


VI 735 Cyg 


21h47m20.6s 


+47d32'00.7" 


3570944 


12 (SL, SH, LH) 


2003-12-17 


V346 Nor 


16h32m32.1s 


-44d55'28.6" 


3570688 


12 (SL, SH, LH) 


2004-02-27 


HL Tau 


04h31m38.4s 


+ 18dl3'57.9" 


3531776 


6 (SL, SH, LH) 


2004-03-04 


XZ Tau AB 


04h31m40.1s 


+18dl3'57.4" 


3531776 


6 (SL, SH, LH) 


2004-03-04 



^The object was observed at three different epoclis. 



Table 3. Overview of prominent spectroscopic features (mostly ices and dust) seen in the spectra presented in 
Figures [H [2], O El El "abs" denotes an absorption feature, "em" an emission feature. For objects where ISO and 
Spitzer data are available, the Spitzer data is of higher quality and thus more reliable. 



OK ^Qr.^- 


Instrument 


ri2U l^icej 


TT pn 


/TTPnmTi'?'! / 


TI2U (^gasj UrLaUrli^.'^j/ 


silicates 


l\r^a\ PAW* 










rL2U i^icej 












o.uo/im 


0. 


/fin iirr, 

oo / o.u /im 


0.0 / o.o/im D.oo/im 


iu.u/im 


LO.Z yfim 


00 Ser 


ISO-SWS 










abs 


abs 


V346 Nor 


ISO-SWS 


abs 




abs 




abs 






Spitzer IRS 






abs 


abs 


abs 


abs 


Z CMa 


ISO-SWS 










abs 




Reipurth 50 


ISO-SWS 


abs 




abs 


abs 


abs 




L1551 IRS 5 


ISO-SWS 


abs 




abs 


abs 


abs 


abs 




Spitzer IRS 






abs 


abs 


abs 


abs 


RNO IB 


ISO-SWS 


abs 








flat 




RNO IB 


Spitzer IRS 


_ 




abs 


abs 


abs 


abs 


RNO IC 


Spitzer IRS 






abs 


abs 


abs 


abs 


VI 735 Cve 


ISO-SWS 


abs 




abs? 




flat 






Spitzer IRS 










abs 


abs? 


Parsamian 21 


Spitzer IRS 










em 


em 


Bran 76 (BBW 76) 


Spitzer IRS 








abs 


em 




V1057 Cyg 


Spitzer IRS 








abs 


em 


abs (gas) 


V1515 Cyg 


Spitzer IRS 








abs 


em 




V1647 Ori 


Spitzer IRS 






abs 


abs 


em 


abs 


FU Ori 


Spitzer IRS 








abs 


em 




V883 Ori'^ 


ESO Timmi2 










abs 




XZ Tau AB 


Spitzer IRS 










em 




HL Tau 


Spitzer IRS 






abs 


abs 


abs 


abs 



Table 3 — Continued 



Object Instrument H2O (ice) H2C0(?)/HC00H(?)/ H2O (gas) CH30H(?)/ silicates CO2 (ice) PAH^ 

H2O (ice) NH+(?) 
3.08/im 5.85/ 6.0 yum 5.8 / 6.8/im 6.85/im 10.0/im 15.2/im 



^Any significant emission at 6.2 or 7.7/im (C-H modes), 8.6, 11.3 or 12.7yum (C-C modes) or at 8.2 /im (origin not clear 



yet). 



Observations from lSchiitz et al.l (120051 ) 



00 
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Table 4. Temperatures and power law indices for the different components of the 
analytical dust model for FU Ori shown in Figure [13l 



Component 


T[K] 


Exponent p 


Inner disk rim 


2246 


-0.84 


Disk surface 


1128 


-0.49 


Disk midplane 


865 


-0.11 



Table 5. Mass fractions of dust species with different grain sizes derived from the dust 

model fit for FU Ori shown in Figure [T3i 



Species 


0.1 fim 


1.5 yUm 


6.0 /im 


Total 


Amorphous olivine 


0.16 


0.20 


0.21 


0.57 


Amorphous pyroxene 


0.00 


0.00 


0.43 


0.43 


Crystalline forsterite 


<0.01 


0.00 


0.00 


<0.01 


Crystalline enstatite 


0.00 


0.00 


0.00 


0.00 


Sihca 


0.00 


0.00 


0.00 


0.00 


Amorphous carbon 


0.00 


0.00 


0.00 


0.00 



Total 



0.16 



0.20 



0.64 
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Table 6. Optical depths and extinction values as derived from the spectra shown in 
Figure [151 The last column list reference values for Ay from the literature. 



Object Tsilicatc Asilicatc [AiHl] ^silicate [mag] Ay [mag] J[Merature ^j^^^g] 



L1551 IRS 5 


1.33±0.03 


9.44 


1.44±0.03 


21.83± 4.70 


19-20^ , 30^ , 


>150^ 


V346 Nor 


0.32±0.01 


9.68 


0.35±0.01 


4.98± 0.86 


6.2'^ 




VI 735 Cyg 


0.29±0.03 


9.52 


0.31±0.03 


4.60± 0.91 


10.0^'^ 




RNO IB 


0.24±0.01 


9.30 


0.27±0.01 


4.15± 1.03 


9.2S 




RNO IC 


0.63±0.01 


9.58 


0.69±0.01 


10.11± 1.87 






HL Tau 


0.47±0.01 


9.61 


0.51±0.01 


7.38± 1.32 




a38j 


Z CMa 


0.31±0.06 


9.48 


0.34±0.20 


5.05± 1.04 


2.8±0.1 


k 


Reipurth 50 


1.68±0.34 


9.68 


1.82±0.20 


26.18± 4.34 


^50^ 





' iSnell fc Schloerbl fll985h 



' iSmith et all fll987h 



ICampbell et all fll988[ ) 



^Graham fc Frogell (11985h . assuming E(B-V)=2.0 



Sato et al.l fll992l ) 



Levreaultl fll988[ ) 



qStaude &: Neckell (Il99ll ). based on optical data and assuming spectral type. 



teergin et all (12005h 



Close et all (Il997f ) 



Men'shchikov et al.l (Il999l ). derived from radiative transfer model. 



1 Cohen &: Kuhil (119791 ). based on optical data and assuming spectral type. 



Casalil (Il99ll ). based on NIR colors, assumed intrinsic color temperatures and strength of 



additional ice absorption features. 
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Table 7. Optical depths of the main ice features at 6.0, 6.85 and 15.2 /xm. 



Object 


76.0 /im 


76.85 nm 


7'l5.2ium 


L1551 IRS 5 


0.51±0.01 


0.47±0.01 


0.54±0.05 


V346 Nor 


0.05±0.01 


0.06±0.01 


0.12±0.02 


V1735 Cyg 


0.05±0.02 


0.09±0.01 


RiO.3 


RNO IB 


0.20±0.01 


0.18±0.01 


0.32±0.04 


RNO IC 


0.35±0.01 


0.26±0.01 


0.41±0.02 


HL Tau 


0.10±0.01 


0.07±0.01 


0.13±0.06 


V1647 Ori 


0.07±0.01 


0.06±0.01 


0.20±0.02 


V1647 Ori 2 


O.lOiO.Ol 


O.OSiO.Ol 


0.19±0.02 


V1647 Ori 3 


0.09±0.01 


0.07±0.01 


_a 



^only low-resolution data 
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Table 8. Upper half: Flux densities at 20 and 30 /im and the slope of the SED between 
15 and 30 /im for 00 Ser observed at 5 different epochs. Lower half: Flux densities at 8, 
20 and 30 fim for V1647 Ori at three different epochs. 



Date Flux density Flux density Flux density Slope [Jy/A] 
at 8 yum [Jy] at 20 fj,m [Jy] at 30 /xm [Jy] 

00 Ser 



1996-04-14 16.5±1 

1996- 10-24 15.2±2 

1997- 03-08 12.8±1 
1997-04-12 11.9±1 
1997-09-22 8.8±1. 



.8 


43.0±2.3 


2.65±0.04 


.1 


41.2±2.6 


2.61±0.05 


.2 


32.5±1.5 


1.98±0.03 


.5 


31.1±1.8 


1.92±0.04 


2 


23.9±1.4 


1.52±0.03 



V1647 Ori 



2004- 10-20 

2005- 03-11 
2005-03-24 



3.95±0.02 
2.49±0.01 
3.04±0.01 



11.56±0.16 

7.21±0.20 

7.82±0.17^ 



16.52±0.51 
10.35±0.24 
10.62±0.05' 



low-resolution data 



V1755Cyg with ISO/SWS 



5 

^ 4 






1 3 




**** * * 


: 2 

D 

^ 1 







6 8 10 

Wovelength [fim] 



12 



14 



L1551JRS5 with ISO/SWS 




6 8 10 
Wovelength [/xm] 

ZCMq with ISO/SWS 



200 
150 
100 
50 



6 8 10 

Wovelength [/im] 



12 14 



RN01B with ISO/SWS 




6 8 10 

Wovelength [/xm] 

ReipurthSO with ISO/SWS 




6 8 10 
Wovelength [/xm] 

V346Nor with ISO/SWS 




6 8 10 
Wovelength [/inn] 



Fig. 1. — ISO-SWS spectra in the 2-15 /xm range for 6 of our targets. The two plots in 
the top panel (V1735 Cyg and RNO IB) suffer from a poor signal-to-noise ratio resulting 
from a combination of relatively low flux levels and a short integration times (Table [1]). All 
sub-structures in the spectra are not believed to be real. 



00Ser_1 with ISO/SWS 




4 6 8 10 12 14 

Wovelength [/im] 



00Ser_3 with ISO/SWS 




4 6 8 10 12 14 

Wavelength [/im] 



00Ser_2 with ISO/SWS 




6 8 10 12 14 

Wavelength [/im] 



00Ser_4 with ISO/SWS 




Fig. 2. — Same as Figure [T] but for the object 00 Ser which was observed at 5 different 
epochs. As shown in Tabled] the integration time of the observations were increased during 
the 5 epochs. However, as the object faded significantly over this period, the signal-to-noise 
ratio could not be improved and it is difficult to identify spectroscopic features in the data 
apart from the silicate absorption profile seen in all data sets. 
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V1735Cyg with ISO/SWS 



RN01B with ISO/SWS 




20 25 30 

Wovelength [/im] 

L1551JRS5 with ISO/SWS 




400 
350 
300 
250 
200 



20 25 30 35 

Wovelength [/im] 

ZCMq with ISO/SWS 



40 



150 



15 



20 



25 30 
Wovelength [/im] 



35 



40 




25 50 
Wovelength [/im] 

ReipurthSO with ISO/SWS 




25 30 
Wovelength [fim] 

V346Nor with ISO/SWS 



50 

3, 40 

"I 30 

■D 

X 20 

_D 

Li_ 

10 









++ + : 












V 







15 



20 25 30 

Wovelength [/im] 



35 



40 



Fig. 3. — Same as Figure [T] but between 15 and 40 /xm. All spikes seen in the data are 
artefacts and not believed to be real, as they showed up either in the up- or in the down- 
scan but not in both. 



00Ser_1 with ISO/SWS 




20 25 30 35 

Wavelength [/im] 

00Ser_3 with ISO/SWS 



40 




20 25 30 35 

Wovelength [/im] 

00Ser_5 with ISO/SWS 



40 




20 25 30 

Wovelength [/irn] 



40 



00Ser_2 with ISO/SWS 




100 



20 25 30 

Wavelength [/im] 

00Ser_4 with ISO/SWS 



40 




15 20 25 30 

Wavelength [fim] 



35 40 



Fig. 4. — Same as Figure [2] but between 15 and 40 /xm. The red, dashed hne illustrates a fit 
to the data between 15 and 30 /i to estimate the change of the slope and flux density levels 
during the decay of 00 Ser (see text section 3.4.). All spikes seen in the data are artefacts 
and not believed to be real, as they showed up either in the up- or in the down-scan but not 
in both. 
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Bron_76 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [fifn] 

VI 647_0ri 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [|tm] 

VI 647_0rr_3 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [fim] 

V1515_Cyg 




0.6 1 , , , , , , , , , J 

5 6 7 8 9 10 11 12 13 14 15 
Wovelength [/»m] 



FU_Ori 



6.0 r 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [fim] 

VI 647_0rL2 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [nm] 

V1057_Cyg 




3 



5 6 7 8 9 10 11 12 13 14 15 
Wovelength [^tm] 

XZ_Tou_AB 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [ixm] 



Fig. 5. — Spitzer low-resolution spectra between 5 and 14 /xm. 
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Parsannian_21 L1551JRS5 




5 6 7 8 9 10 1 1 12 13 14 15 5 6 7 8 9 10 1 1 12 13 14 15 

Wovelength [/xm] Wovelength [fim] 



V1735_Cyg RN0_1C 




5 6 7 8 9 10 1 1 12 13 14 15 5 6 7 8 9 10 1 1 12 13 14 15 

Wovelength [/im] Wovelength [/im] 



RN0_1B V346_Nor 




5 6 7 8 9 10 1 1 12 13 14 15 5 6 7 8 9 10 1 1 12 13 14 15 

Wavelength [/nm] Wavelength [/um] 



HL_Tou 

lOf ' ' ' ' ' 




5 6 7 8 9 10 11 12 13 14 15 
Wovelength [/xm] 



6. — Same as Figured 
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Brfln_7B 




u m ia 20 23 24 za U 34 




1^ ri ID 20 2B 34 2E £B » J2 34 
IhiyBlBiiglh Ij^j 

Vl&4?_0fi_a 




14 ri IS 20 22 :4 S IS X 32 34 



vnSi5jCy 




14 T| IB JO 22 3* 38 IS 30 K 34 




M II ia 20 33 2>t 2fi 29 30 32 3* 





14 1S i« 2D 2? 2S a 10 3S 3* 




14 11 ie 20 Z7 24 » 29 3S ^ 3« 



Fig. 7. — Same objects as in Figure O but now showing the high-resolution part of the 
spectrum between 14 and 34 /xm. The spectrum of V1515 Cyg shows a prominent emission 
hne longward of 33 /im which is discussed in section 3.5. Between 31 and 32 fim all objects 
where high-resolution data is present show artefacts from the data redcution process. 
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Ptlf«ffliOrt_2l L155t_:R55 





I 1 

r 



U 11 ZD 23 34 SB U M 




14 re IB » 22 3* » 3a JO 32 3^ 



Hl_lau 




M IB m 3S 23 2* 2i 2a 30 32 » 



RM0-1C 




M IS ^fl 30 ;^ ZS JQ JZ 9<. 




14 19 ie 29 22 74 S» 2S a J2 3« 



Fig. 8. — Same objects as in Figure [6l but now showing the high-resolution part of the 
spectrum between 14 and 34 /xm. Between 31 and 32 fim all objects where high-resolution 
data is present show artefacts from the data redcution process. 
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L1551 IRS 5 



V346 Nor 




7 8 9 10 11 12 
Wavelength [/xm] 

RNO IB 




7 8 9 10 11 12 13 14 15 
Wovelength [/xm] 




8 9 10 11 12 13 14 15 
Wavelength [fJ.m] 



VI 735 Cyg 




8 9 10 11 12 13 14 15 
Wavelength [jitm] 



Fig. 9. — Comparison between ISO/SWS and Spitzer/IRS data for objects observed with 
both sateUites. The ISO data (red, sohd hnes) are significantly more noisy. In our sam- 
ple only data for objects with high flux levels or long integration times are suitable for 
quantitative analyses. 
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Fig. 10. — Continuum subtracted silicate emission profiles (black, solid lines). Each spectrum 
is shown twice, i.e., with and without being corrected for interstel lar extinction as expla ined 
in the text. T he assumed values for Ay are: 2.2 mag ( Bran 76: iReipurth et aL 2002 ). 2.4 
mag (FU Or i; Iskinner et allbooeh . 3.2 mag fV1515 Cve: EeTb igJIlOTTh, 3.2 m ag (V1057 Cyg; 



Herbidll977h . 11.0 mag (V1647 O ri: iMuzeroUe et a 



Iboosh . 1.4 mag (XZ Tau; IWhite fc GhezI 
1992h . We used the extinction law from 



20011). 4 ma g (Parsamian 21; IStaude fc Neckell 
Mathid (jl990l ) with Rv=3.1. For comparison the red, dashed line shows t he typical ISM 



silicate feature observed toward the galactic center from lKemper et al.l (120041 ). As mentioned 
in the text V1647 Ori was observed at three different epochs. 
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8 9 10 11 12 13 8 9 10 11 12 15 

Wavelength [fj-m] Wovelength [/xm] 



Fig. 11. — Dereddened and normalized 10 /xm dust emission features of the objects shown in 
Figure [101 For a better comparison the dereddened spectra are shown. 
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K4 - 



K2 - 



1,0 - 



0,8 - 



0,6 - 



8 



1 2g 5 
^3 



9 
X 



1: FU On 
2: Bron 76 
3: VI 51 5 Cyg 
4: VI 057 Cyg 
-5:^1647 Ori 
6:"vi'B47 Qr\^2 
7: VI 647 Orr3"". 
8: XZ Tou 
9: Porsomian 21 



1,0 1,5 2,0 2,5 

Peak over continuum [nnax(F„orm)] 



Fig. 12. — The flux ratio between 11.3 and 9.8/im plotted against the normalized peak 
over the continuum. The dashed line indicates the region where typ i cally TTauri stars 
(jPrzygodda et al.ll2003l ) and Herbig Ae/Be stars (Ivan Boekel et al.ll2003l . l2005l ) are located. 
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Fig. 13. — Results from a dust model fit to the 7-17 /im region of the spectrum of FU 
Ori using the model described in the text. The observed spectrum is shown in the black, 
solid line, the blue, dotted hne illustrates the contribution from the fitted continuum, and the 
green, dashed hne denotes the computed emission feature. The sum of the fitted components 
is shown in the red, dash-dotted line, matching the observed spectrum. The mass fractions 
of the fitted dust species are given in Table [51 The spectrum is fitted solely with amorphous 
grains, part of which have grown significantly. 
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4 6 8 10 12 14 16 
Wavelength [/xm] 

Fig. 14. — Dereddened Spitzer MIR spectrum of Parsamian 21 (assuming Ay = 4.0 mag) 
illustrating the strong PAH emission features not seen in any other FUor. Note also the 
steeply rising continuum and the underlying silicate dust emission. The vertical lines indicate 
the positions of typical PAH emission bands: strong emission is detected at 6.3, 8.2 and 
11.3/im (red, dashed lines), no emission is seen at 7.7, 8.6 and 12.7/im. 
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Fig. 15. — Observed optical depths of the 10 /im sihcate absorption feature (black, solid 
line). The dotted line shows the assumed continuum. The vertical red, dash-dotted lines 
indicate the position of the maximum optical depth of the second order polynomial fitted to 
the data between 8.8 and 10.2 (also plotted in red, dash-dotted lines). 



-57- 




8 9 10 11 12 15 
Wavelength [fxm] 



Fig. 16. — Optical depths for different dust grain populations on an arbi trary scale. The 



black , solid line shows the profile for 0.1 fim sized amorphous olivine grains (IDorschner et al. 
19951). the red, da shed line shows the ISM feature toward the galactic plane from 



Kemper et al.l (12004). the green, dash-dotted line is the profile based on the astronomi- 



cal silicates from IWeingartner fc Draind (1200 ll ) andlDraind (20031). an d the blue, dotted line 
is for 0.1 fim sized amorphous pyroxene grains (jPorschner et al.lll995l ). 
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Fig. 17. — Same observed spectra as in Figure fT5l (black, solid lines), but now overplotted 
with reference spectra with different dust compositions (same color and line style code as in 
Figure [16]) . 
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Fig. 18. — Compari son of the 15.2 /xm C02-ice feature between the heavily e xtincted back- 



groun d object CK2 (IKnez et aLll2005l ) and the embedded young source HH46 (IBoogert et al. 



20041 ). Due to higher temperatures in the vicinity of the embedded protostar crystaUization 



occurs and and a double- peaked sub-structure appears characteristic of the pure CO2 matrix 
( lEhrenfreund et al.lll998l ). 
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Fig. 19. — The 15.2 /iin C02-ice feature observed toward some of our objects (black, solid 
lines). The spectra of HL Tau, RNO IB and V1735 Cyg were smoothed by a factor of 20 to 
increase the signal to noise. Overplotted is either the scaled spectrum of HH46 or CK2 (red, 
dashed-dotted line; see Figure [18]) depending on which profile fits better to the FUor data. 
Due to the remaining high noise level in the spectrum of V1735 Cyg no clear comparison to 
either reference object was possible. The absorption feature of V1647 Ori observed at the 
third epoch is not shown as these are only low-resolution data, while the data presented here 
(first and second epoch) were taken with the high-resolution spectrograph. 
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Fig. 20. — The 6.0 and 6.85 /im ice bands observed toward some of our objects (black lines). 
Overplotted are the polynomial fits of fourth order (red, dashed-dotted line) to determine 
the optical depths of the absorption bands. 
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Fig. 21. — Zoom into emission lines and absorption lines detected in the high- resolution part 
of the spectra of some objects (see text). 
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Fig. 22. — Correlation between the optical depths of ice features at 15.2 and 6.0/im (left), 
ices and silicates at 15.2 and 9.8/xm (middle), and ice features at 6.8 and 6.0/im (right). 
V1647 Ori is not shown in the middle plot as its silicate feature is seen in emission and not 
in absorption (see text). 
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Fig. 23. — Spectral index n of the silicate emission objects evaluated between 13 and 25 /xm 
vs. the spectral index between 6 and 13 /im. n is computed as n = (ilog(AFx)/(ilog(A). 
AH objects populate t he same area in the plot as the Class II objects in Taurus presented 
by iFurlan et al.l (120061 ) (dashed area). Due to the steep rise in its SED, Parsamian 21 lies 
slightly off from the other objects. 
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Fig. 24. — Simple graphical representation explaining the two categories of FUors (not to 
scale). Category 1 objects (left) show silicate absorption and are younger than than Category 
2 objects (right) which show silicate emission. While Category 1 objects are still more deeply 
embedded in their envelopes and appear to be Class I sources, Category 2 FUors are similar 
to Class II sources. The spectrum of V1647 Ori can be explained assuming that extinction 
in the line of sight to this object creates the ice absorption features and reduces also the 
strength of the observed silicate emission. 



